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SEMI EMP I RI CAL METHODS FOR COMPUTING TURBULENT FLOWS 
I. A. Belov, I. P. Ginzburg 


Determination of turbulent friction and heat exchange /159 * 

in a wall region is based on two semiempir ical theories: the 

Prandtl-Karman mixing path theory and a theory using equations 
for turbulent pulsation energy with additional Kolmogorov-Rotta 
hypotheses. 


1. Prandtl-Karman Theory and Effective Viscosity 
Hypotheses. Existing models for Prandtl-Karman ' s 
semiempir ical theory of turbulence applied to boundary layer 
problems are distinguished by various assumptions regarding the 
size of the mixing path and the presence of a laminar 
sublayer. They underlie hypotheses that, when flow at the 
boundary layer of an axisymmetric wing or body is in question, 
turbulent flow may be calculated with Prandtl's equation 


t = P l 2 


d*>x 

dy 


dVr 

~dy~ ’ 


where 1 is the mixing path. 


( 1 . 1 ) 


According to the two-layer system, it is assumed that the 
boundary layer may be broken down into a laminar sublayer in 
which laws of viscous friction and heat transfer are valid and 
a turbulent layer in which Prandtl's equation is valid. The 
thickness of the laminar sublayer 6 -^ is determined from the 
condition of a break in the derivatives of velocity at the 
boundary of the laminar sublayer. 


dy 


ly=8. -0 


. v- dV x 

K '~dT 


ly-*/+o 


( 1 . 2 ) 


*Numbers in the margin indicate pagination in the foreign text. 
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where K-^ is a certain universal constant. If one assumes 
that 1 = Ky near the laminar sublayer, the following formula 
for calculating results: 



(1.3) 


This model was consolidated for a compressible fluid and a gas 

mixture. In addition to the two universal constants K and 

K^, we must introduce the concept of mixing paths for 

temperatures and for individual components in the mixture or, 

T T 

respectively, Pr = 1/1 T and Sm^ = 1/1^, where 1^ 
and 1^ are mixing paths for temperatures and for 
concentrations of the i-th component in a mixture. 


As numerical analysis has shown, in calculating friction /160 
and heat transfer the hypothesis 1 = Ky gives satisfactory 
results which correlate with experiments. The universality of 
constants K and produces the limit relationships derived 
by Kutateladze and Leont'yev [1]. To determine a velocity 
profile in a boundary layer which correlates with experiments, 
one must make more general assumptions about the mixing path, 
taking into account the fact that it has a constant value near 
the edge of the boundary layer. Here we might point out 
various functions 1/(5 = f(y/<T) proposed by several authors. 
Spalding's standardized formula [6] was used at Leningrad State 
University (LGU) to calculate 1/6 in numerical analyses of a 
turbulent boundary layer. 




(1.4) 


where is a certain constant close to 0.1; K = 0.4. 


In a three-layer system (the Klauzer-Mellor hypothesis 
[7-8]) the boundary layer is divided into three regions: a 

laminar sublayer and the inner and outer regions of the 
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boundary layer. It is presumed that: 


y<*i, 


—nr -Ss.. 


< 


dOjf 

^IT’ 

So Uh* 


5,: 


K t 


K ^ v ^ K v 

is the boundary layer interior; 


K /x/p . 


:p /^2 y 2 J^r 



Ui* 


T = p o pf7o* 


0v r 

1y~ 


is the boundary layer exterior, where 6* is displacement 
thickness; |3 q is a certain constant (according to Klauzer, it 
equals 0.016-0.018; in calculations done at LGU, it was set at 
0.02); U is the velocity at the outer edge of the boundary 
layer. 


There are also studies where the boundary layer is not 
divided into a series of regions, but a so-called effective 
viscosity is introduced, and it is assumed that t=i i e f(d/dy)v x , 
whereby Pef = P+PT, , where, according to Prandtl p T = p l 2 (d/dy)v x . 


Mixing path 1 is determined differently for different 

regions of change in y. To account for the fact that the 

3/2 

mixing path near a wall approaches zero at y , van Driest' s 
equation [9] is used to calculate it in the boundary layer 
interior. 

is =Ky[l — exp(— '^r)/ r ~)]. T^X’ (1.5) 

where A is a certain constant (according to van Driest, A = 

26). According to Spalding [6 ] , in the boundary layer exterior 

i = K i X =0,1; AT=0,4. (1.6) 
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To account for the effect of blowing and a pressure gradient /161 
for an incompressible fluid, Sebesi [10] assumes the following 
for the boundary layer interior 

/=;<>[i_exp (1,7) 


whereby ~f~ w is determined from an approximated solution to the 
equation 




tx 


dp i <*t 
dx 1 dy ’ 


V„ 




( 1 . 8 ) 


from which it follows that 



In the case of a compressible fluid, (1.8) gives 


(1.9) 




dy + t w 


exp 



In the absence of blowing. 



(1.9*) 


and equation (1.7) takes the form 


/ = /0/[l-exp(--^-|A 


(1.7') 


With regard for blowing, 


^^ = fe^(-> + e»p-,) + ex P ^- 8< r, (1 . 10, 


wherein Sebesi assumes 5 /=^--^r= H,8 where v*=V \JP- 

In Sharov and Lapin's work [11], the effect of blowing and a 
pressure gradient was taken into account on the assumption that 
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2 


( 1 . 11 ) 


0 ,= 

Ki v f dp 
K v * [ dx 


K\ v __ K x * ]/ %_ 

^ Tj,/p K ' "‘I 


-1— f— 1 + exp 

pl'a, \ K 


(*• 


»/) + «P- £ t ! -8|' 


This equation is solved numerically. Note that studies by 

Sebesi, Lapin, and Sharov divide the boundary layer into 

interior and exterior parts. In the interior, 

where 1 is calculated according to (1.7) with regard for 

(1.10); in the exterior, jXj, is determined from Klauzer's 

equation 


Velocity profiles are joined on the basis of the equality 
of turbulent viscosity and of velocity continuity coefficients 
and velocity gradients. 


To evaluate the possibility of using particular models of 
the semiempir ical theory of turbulence, LGU numerically solved 
equations for movement of an incompressible fluid in a boundary 
layer in the presence of blowing and a pressure gradient 
(certain calculation results have been published in articles by 
G. V. Kocheryzhenkov, S. K. Matveyev, and V. S. Ivanov [2-5]). 

These calculations showed that discrepancies between results 
for one system and for the other in functions for friction 
coefficients and for K = &*/&** from x, where x is 

the distance along a plate, are within limits of calculating /162 

accuracy. Figures 1-3* present certain results of the 
comparison. Note that numerical analysis is naturally simpler 
in the two-layer system, where there is no need to formulate 
clarifying hypotheses regarding the mixing path. 

Along with numerical solution of equations in partial 
derivatives, approximated integral methods of calculating 
friction and heat transfer in the case of a turbulent boundary 
layer are used. These methods are based on use of integral 


*Figures are taken from [2-4] . 
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relationships from boundary layer theory and an approximated 

concept of friction stress in the form of polynomials of y /6 

and v /U such as 
x 

t/^==l +-4-"} — ( ! +4") y ‘ 2 ' (1.12) 

(1.13) 

where ?I_=y/ 8 , \-=v x jU. 



Fig. 1 Change in surface friction 
factor . 

1, 2 - Numerical solution 
according to Mellor's hypothesis 
and the two-layer system; 3, 4 - 
Approximated solution using the 
same models. 


Figure 2. Change in H = 



Figure 3. Surface friction factor in the turbulent 
boundary layer of a plate during blowing: 
a - v ^/\J - 0.00386; b - v^U = 0.0095; 1 - Two-layer 
system; 2 - Effective viscosity system; 3 - Sebesi's data; 
4 - Simpson, Kays, Moffat experiment. 
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From conditions taken from motion equations, assuming that 
they are valid also near a wall, we find that 


dp U\y w i ?w v ki U 

dx 

w 

dp JS , pt/ TO S 

dX t w ‘ Ra; 


(1.14) 


Comparison of numerical analyses with approximated solution (if /163 
A in (1.12) is calculated from (1.14) when a pressure gradient 
and blowing are present) showed considerable discrepancy among 
results. This is because the effect that friction distribution 
over the surface has on the equation for A is not considered. 

In studies by G. V. Kocheryzhenkov and S. K. Matveyev [2, 12], 

A was calculated from the condition A — — t/t„, I 
In this case. 


dp 

dx 


+ 


_L dv * ( v * W* 

v* ' dx \ U j p 




Ki 


exp 


[*•-*%*) 


(1.15) 


where • During calculation, it was assumed 

that ^ = 0.6. With regard for A in terms of (1.15), 
approximated calculations coincided with numerical solutions. 
The approximated solutions produced on the assumption of the 
law for T/T w (1.13) and for B (1.14) showed a rather good 
correlation with the numerical solution. 



Fig. 4 Profiles of longitudinal 

root-mean-square velocity pulsation 

in the boundary layer with external 

flow turbulence £ = 4.6%. 

a 

1 - Calculated using van Driest' s 
ratio for mixing flow; 2 - Same as 
per Spal ding ; 3-5 - Experiment; 
i — Jv^/U = 0.049; 

Re x = u x/v = °* 184 x 1q5; 

4 - 0.072; 10; 5 - 0.083; 1.53 

x 10 5 . 
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In summary, one might note that this comparison of the 
results of calculations for different systems shows that 
calculating friction and heat transmission requires numerical 
computation of boundary layer equations and determination of 
turbulent friction for the two-layer system. Approximation 
methods should be used in evaluation, assuming x/x w =f( I) in 

accordance with (1.13) and (1.14) . To find thermal flows, one 
can use well-known relationships between heat content and 
velocity as applied to gas mixtures [13] . 

2. Semiempir ical Theory of Turbulence Involving Use of 
the Equation for Turbulent Pulsation Energy. In solving the 
problem of impingement of flows with external turbulence 
(natural or imposed) upon a barrier, it is essential to account 

for external flow turbulence when friction and heat exchange- 

with the barrier are being analyzed. It is more convenient to 
solve this problem using an equation for turbulent pulsation 
energy, since, in addition to friction and heat exchange, it 
permits determination of turbulent pulsation energy with regard 
for external flow turbulence. Note that the hypothesis on 
mixing flow is ' unacceptable here, since it gives zero values 
for pulsations at the outer edge of the boundary layer. An 
example of a boundary layer with external turbulence is the 
boundary layer at a plate immersed in a jet flow. The results 
are given below (fig. 4) [14] for experimental study of such a 

boundary layer when a plate is immersed in a subsonic jet 
within the limits of the initial segment of the jet 
(gradient-free flow) with turbulence intensity at the nozzle 
exit section 

/ u] =A, 6 % 

As follows from data in this figure on distribution of 
longitudinal (along the su rfac e of the barrier) root-mean- /164 

square velocity pulsation^' ^ in the barrier's boundary 
layer, velocity pulsations in the jet penetrate deep into the 
boundary layer and create conditions for acceleration of the 
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transition from laminar flow (near the nozzle exit section) to 
turbulent. Maximum velocity pulsations occur in close 
proximity to the wall: in the' laminar flow area, at a distance 

of about y/4 = 0.2-0. 3, where & is the thickness of the 

boundary layer; in the turbulent flow region, at a distance y/cf 
<0.05. The f igure also presents results of computing velocity 
pulsation^' 2 = ld/5y(v ) , where v /v* = 5.5+5.75 lg 
y v*/v; v* 1 is the mixing path as per van Driest 

(curve 1) and as per Spalding (curve 2) . As these data show, 
the results of computation correlate well with the experiment 
only very close to the wall. This fact raises doubts about the 
existence of a linear relationship between velocity pulsations 
and the gradient of velocity averaged over time in the 
turbulent boundary layer's turbulent region. 


The semiempir ical turbulence theory based on use of energy 
equations for turbulent pulsations makes it possible to define 
all flow parameters in greater detail. However, to close the 
complete system of equations requires we must introduce several 
additional coefficients which are not always known and not 
always universal. Let us evaluate these coefficients. The 
energy equation for -turbulent pulsations for an incompressible 
fluid takes the form 



C2.1) 


where v^ is the constituent in direction x^ for motion 
velocity averaged over time; v* . is the pulsation of v.; 

1 rjT ^ ^ 

£=Y Zj Vn • . It is usually assumed that 

n 




n 


( 2 . 2 ) 



.1* 


where o' is Prandtl's turbulent number for turbulent pulsation 
energy; is turbulent viscosity; c D is a constant; 1 is 
the scale of the turbulent pulsations; 5 — is Kroneker's 
tensor. 


With regard for (2.2), equation (2.1) is rewritten as 



(2.3) 


In accordance with Prandtl-Kolmoqorov's equation, turbulent /.'L65 

1/2 

viscosity is expressed at v T = c v E ' 1. One of the most 

important relationships in (2.3) is that between the kinetic 

energy of turbulent pulsations and energy dissipation. 

According to Kolomogorov's hypothesis on local isotropy, c D in 

(2.3) can be considered as a universal constant when the subject 

flow is characterized by very high Reynolds numbers for 

1/2 

turbulent pulsations, Re^ = E 1/v. Given low and moderate 
Re T values, the lower the Re T value, the more c D deviates 
from its limit value. Rotta (15] proposed the following 
approximated function for c^s 

(2.4) 

From (2.4) it follows that the limit value for c D is 0.202 as 
Re T approaches infinity. According to Wolfstein's data [16] 


From data in [17] , 

c D = Re T /(/?), 


= _ — y A n =0,263. (2.5) 

1 — exp (~A d Re T ) > D 
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where 


/(R) = r*^.'»r(-T. -i-) + fi ! ' 3 «p(- #)]■’• 
c = 1.5; x s= 4/3; a = 7.5; f" is a gamma function; 

R = CD m l 5:4 /{^ 4 E 3i8 ), ( 0,2 </?< 1000 ). 


It follows from the last function that, when Re T approaches 
infinity, c D approaches 0.15 instead of 0.202. At low Re T 
values, c D increases more quickly than is indicated in (2.4). 


To find the connection between c and c in equations 
for turbulent viscosity and turbulent energy dissipation, let us 
study the particular case of equation (2.3), which described 
flow of a fluid in a state of local equilibrium [18] . For this 
flow, equation (2.3) is rewritten as 


But 


from which 


cJ ;’° l (&y = c D E* 


% 


v r v 


— —cE v ‘ l l 

^ur~ ' 

dy 


p£ 


g - — ( c £) c v ) 


1/2 


According to data in [18] 

t/p£'.2=s0,32 

for a uniform transverse flow and for flow in tubes near a wall, 
but still rather remote from the viscous sublayer (one may 
approximately state that these flows are in a state of local 
equilibrium) . Consequently, c c v Or 0.1. Assuming, as did 
Ng and Spalding [18], c Q = 0.1, we obtain c = 1.0. Taking 
c n = 0.416, as did Wolfstein [16], we obtain c = 0.24 
(according to Wolfstein's data, c v = 0.22). If c v does 
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not equal 1, assuming c' to be constant for the entire range 
of change in Re T , Wolfstein recommends the following equation 
[16] s 

r, = a v [l — exp(— ARe T )], A-, = 0,016. (2.6) 

Prandtl's turbulent number incorporated into (2.3) varies /1 6 6 
within 1.53-2.5 [15, 16, 18]; the nature of the change in scale 
1 in (2.3) must be determined in each specific case. 


As an example of the use of the turbulent pulsation energy 
equation, let us calculate the impingement of a flat turbulent 
jet to a plate set along the normal to the jet flow. This flow 
is characterized by the existence of a stagnation point on the 
plate. Directions along the plate (x) and along the normal to 
it (y) for the area around the stagnation point are 
equivalent. Therefore, study of this flow requires complete 
Reynolds equations. These are written in the form of a vortex 
transfer equation as 


v 


dQ 
* dx. 



(2.7) 


To equations (2.7) and (2.3), we will add an equation for 
continuity 


dv x . dv y _ n 
dx ■~ r ay 


( 2 . 8 ) 


With regard for the relationship of the components of the 

Reynolds stress tensor (2.2) and the Prandtl-Kolmogorov 

equation, the system of equations (2.3), (2.7), and (2.8) is a 

system of three equations for determining four unknown 

functions: v , v , E, and 1. Therefore, one more 

x y 

relationship is required to close this system. For an ordinary 
turbulent boundary layer on the surface of a solid body, this 
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relationship is the familiar mixing path equation 1 = Ky, K = 
0.4 or any other of the mixing path equations presented above. 
If there are pulsations at the outer edge of the boundary 
layer, this approximation may be used as the first 
approximation. 


The velocity profile in the jet remote from the barrier 
(unperturbed flow) in a direction normal to the barrier is 
represented as Schlichting ' s profile. Then, assuming that the 
velocity constituent normal to the barrier in the area of 
interaction with the barrier varies linearly with distance from 
the barrier, we obtain the following directly at the barrier 
[19] if 0 < x i r: 

UjVv = b+ V^l'-y^ + T^j- 

V/Vv = — ( 1 — C 3/2 ) 8 y Vih, - (2 * 9) 

where 4 = (x-b)/6 c ; r is the free jet's boundary; b is the 
boundary of constant velocity nucleus; & is mixing zone 
thickness; Y is the velocity gradient at the stagnation point; 
v is the kinematic viscosity factor. Equation (2.9) is taken 
as the boundary condition at the outer edge of the viscous 
mixing of flow and plate yV y/v-*-°°. . Let us formulate 

boundary conditions for the outer edge of the viscous mixing 
layer for turbulent pulsation energy. Flow far from the 
barrier beyond the limits of the boundary layer' may be 
considered with a sufficiently high degree of accuracy to be 
identical to flow in a jet stream for which [20], by analogy 
with an axisymmetric jet. We assume: 


Then 


= +' y,2 *] c - 


( 2 . 10 ) 


VlA, - V t»i - V v* = X ( v ~ , 


/167 
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where 


V^v'l c , Vv'yc, V v'h 

are the velocity constituents of pulsations in a free jet in 
the direction of the axes of a Cartesian coordinate system x, 
y f z, respectively; Y is the distance from the nozzle exit 
section to the barrier; Y*o is the distance from the barrier to 
the cross section in which the effect of the barrier on the 
flow is negligible (as per [12] ) f Y,*, ~ d , where d is the 
diameter of the nozzle's outlet section; V is the free jet 
constituent normal to the barrier (it is described by 
Schlichting ' s profile); x is a constant (as per [20], x = 
0.0256). 


Equation (2.10) assumes a zero value for pulsation energy 
at the axis and outer edge of the jet, which corresponds to 
reality for small distances between nozzle exit section and 
barrier Y. As Y increases, the equations for E c suggested in 
[22] are preferred 


0,22 v\ 

E *~ (Y - y<o)ld a ‘ 


( 2 . 11 ) 


Upon impingement of a plate, pulsation energy changes in the 
area where the jet interacts with the plate (outside the limits 
of the flow/piate mixing layer) . Therefore, (2.10) or (2.11) 
may be used only in the first approximation as the boundary 
condition for pulsation energy at the outer edge of the 
boundary layer. 


Boundary conditions on the plate's surface follow from the 
equation for fluid adhering to that plate at y = 0: 

v x -v y = E=0, y=z 0. (2.12) 


We seek the solution for v x , v^, 
law of change in 1 in the form 


E given the well-known 
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(2.13) 


<x> 3 (m — 1) 

vJVv = rV^'2 i (-l) m+1 C 


m=*l 


3m— 3 


VV? = - (- ')‘ +1 • ■' 


m =* I 


Taking 


so- 


where 


P (C) = (l — C 3 ' 2 ) 2 , 1 c — * 


the solution for E takes the form 


£-• fr°< C) 


s<-» 


3m —3 

m+l c 2, c 


m— 1 


-r 


(2.14) 


where f , S are functions of J subject to definition from 
m m 

equations for vortex transfer and turbulent pulsation energy 


“ — T*( K *”* V{Y~ YMi ' 
dG,=(;>-t«)*/(i 

The last function is present in the form 

. G (C) = 1 -f A 1 t 3/2 + A S C 3 + AJ* 12 + AJ + A t C 15/J -f .4/. 9 + A ? ; 21 ' 2 + ... > 

where = -0.4; A 2 - -0.78; A^ = -0.392; A^ = -0.0395; 
A c = 0.2231, A = 0.3033; A-, = 0.0061, etc. 

DO/ 


2,35 V a 
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After substituting (2.13) and (2.14) into (2.3) and (2.7), 

equation coefficients with identical degrees of , we obtain a 

system 2m of ordinary fourth-degree differential equations. 

Boundary conditions for functions f , S follow from 

conditions (2.9) -(2.12) when y>f y/v = 0 and y*>f = °°. The 

system of equations obtained in this way was solved numerically 

using the Runge-Kutt method on a BESM-6 computer with an 
-5 

accuracy of 10 by sequential approximation. The number of 
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Fig. 5. Profiles of longitudinal 
average velocity in the boundary 
layer of a plate when a subsonic 
flow impinges upon it along the 
normal: 1 - £ = x/r = 0.1; 0.2; 2 

- 0.3; 3 - 0.4; 4 - 0.5; 5 - 0.6; 
0.8; 6 - 0.7. 

Fig. 6. Distribution of friction 
over the surface of a barrier when 
a subsonic jet impinges upon it 
along the normal; 1 - Current 
analysis; 2 - data from [109] with 
x = 0 . 


approximations is determined from the condition 

|B m '(oo, Oj<s=10-V 

where 0 = f or S ; t, t+1 is the number of 

m m m 

approximations; £ is the given accuracy with which the boundary 

conditions are fulfilled at y^y/v -> CK> . Additional boundary 

conditions required to begin calculation were found with 

Newton's method to satisfy boundary conditions for f , S at 
mm 

y'fy/v-^with given accuracy £.. 

Figures 5 and 6 present the results of calculation with the 

-5 

following raw data: V a = 20 m/sec; v = 1.5*10 

m 2 /sec; d 3 = 10 mm; Y/d a = 8; V = V /d, [21]. The 

a a a a 

following approximations were used in calculation: 1 = Ky, K = 

0.4; CT = 1.7; c D = 0.15; c y = 0.66. Equation (2.10) was 
used as the boundary condition for pulsation energy at the 
outer edge of the boundary layer. 

Review of these graphs shows that velocity profiles (fig. 
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5) in various cross sections along the plate are typical 
profiles for the turbulent boundary layer when the pressure 
gradient for the exterior flow is negative. The boundary 
layer's thickness initially drops as it moves away from the 
stagnation point, then increases. Maximum thickness 
corresponds to the cross section with maximum pulsation energy 
at the outer edge of the boundary layer. The solution shows /I 6 9 
that exterior jet flow turbulence has an effect on flow in the 
boundary layer throughout the subject area, including the area 
around the stagnation point. Friction at the plate (fig. 6) 
quickly increases due to the negative pressure gradient in the 
exterior flow, then decreases near the jet boundary due to 
dissipation effects. One might note the fact that friction at 
the plate for this type of turbulent flow-plate interaction 
increases the corresponding figure [19] , obtained if the effect 
of the turbulent jet on flow in the boundary layer around the 
stagnation point of the plate, is disregarded by a factor of 
1 . 8 . 

List of symbols used: x, y, x -- axes of the Cartesian 

coordinate system; 1, 1 T , 1^ -- mixing path, mixing path 

for temperature, for the i-th component of the mixture; S, 

-- boundary layer and laminar sublayer thickness; b --the 

boundary of the core of constant velocities for a free jet; r 

-- the free jet's boundary; Y --distance from the nozzle exit 

section to the barrier; Y -- distance from the barrier to the 

cross section at which the effect the barrier has on flow is 

negligible; v , v , v -- components of velocity toward 
x y z. _____ > . 

axes x, y, and z averaged over time; j/?Jy 2 , .}/"££* 

root-mean-square pulsations of velocity components v x , v , 

v ; U, V -- velocity components toward x and y on the outer 
z 

edge of the boundary layer; v^ --dynamic velocity; £-- 
turbulence intensity; Y -- .velocity gradient at the stagnation 
point of a flat barrier set along the normal to the jet; p , p 
-- density, pressure; T -- friction stress; 

Ju, JLL Q £ - coefficient for dynamic viscosity, turbulent 
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viscosity, effective viscosity; E -- turbulent pulsation energy; 7j = y/8 

%~v x !V\ C = (x — 6)/3 c ; K, Ki — universal constants in (1.3); 

P.= (6) y .y Po -- constant in Klauzer's equation; A, B -- 

coefficients in the breakdown of T/T w (1.14); c D , c -- 

constants in (2.1); x -- constant in (2.10); f , S 

m m 

functions of y defined from (2.7) and (2.3); t -- approximation 
T T 

number; Pr , Sm^ -- Prandtl's, Schmidt's turbulent 

number; (f -- Prandtl's turbulent number for turbulent pulsation 

energy; Re T -- Reynolds number for turbulent pulsations; 

<^j -- Kroneker's tensor. 

Subscripts: i, j, n = 1, 2, 3; m = 0, 1, ..., ; a -- 

parameters at the nozzle's exit section; c -- parameters in a 
free jet. 
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